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We  developed  a  technique  using  nonlinear  correlation  of  photoluminescence  (PL)  to  characterize 
midwave  infrared  lasers  by  extracting  the  density  and  temperature  dependence  of  the  carrier  lifetime 
and  its  exact  branching  into  radiative  and  nonradiative  processes.  This  was  accomplished,  without 
time  resolving  the  PL  recovery ,  through  mathematical  optimization.  We  extracted  this  information 
by  using  a  laser  source  that  can  be  operated  in  both  continuous-wave  and  short  pulse  modes. 
Through  fitting  of  the  PL  signal  and  its  nonlinear  correlation  for  both  laser  modes  of  operation,  the 
carrier  lifetime  as  a  function  of  density  is  extracted.  As  a  proof  of  principle,  we  investigated  a 
midinfrared  Sb  based  laser  and  showed  that  the  radiative  branching  ratio  drops  from  —54%  at  80  K 
to  about  3%  at  room  temperature,  resulting  from  an  order  of  magnitude  increase  in  the  nonradiative 
rate  coupled  with  a  factor  of  2  reduction  in  the  radiative  rate.  We  believe  that  this  is  a  very  generic 
approach  and  can  be  extended  to  various  luminescing  material  systems.  ©  2008  American  Institute 
of  Physics.  [DOl:  10.1063/1.2894310] 


I.  INTRODUCTION 

Midwave  infrared  (MWIR)  lasers  and  detectors  operat¬ 
ing  in  the  range  of  2-5  pt m  are  in  great  need  for  wide  range 
of  applications  such  as  remote  sensing,  LADAR,  detection  of 
chemical  warfare  agents,  intelligence,  surveillance  and  re¬ 
connaissance,  missile  tracking  and  infrared  countermeasures. 
The  most  promising  technology  for  fabricating  MWIR  lasers 
employs  antimonide  based  materials.  Electronic  and  opto¬ 
electronic  devices  based  on  GaSb  substrates  have  been  a  sub¬ 
ject  of  active  research  in  the  past  few  years. 1-4  In  the  past,  Sb 
based  optically  pumped  lasers  have  been  reported  on  GaSb 
substrates  using  a  metamorphic  buffer  grown  with  digital 
alloys.5  These  are  “As-free”  lasers  and  are  based  on  a  type-I 
alignment  in  the  InGaSb/InAlSb  system.  Using  a  digital  al¬ 
loy  buffer  layer  terminating  with  Al  ^In^vSb  as  the  “sub¬ 
strate,”  a  multiquantum  well  active  region  using  AlGalnSb 
barriers  and  Ga6In4Sb  wells  demonstrates  both  photolumi¬ 
nescence  and  intense  cw  laser  emission  at  2.5  /rm  under  op¬ 
tical  pumping  by  a  980  nm  pump  laser  at  room  temperature. 

Accurate  characterization  of  the  recombination  pro¬ 
cesses  in  these  structures  is  indispensable  for  guiding  and 
optimizing  the  growth  process  toward  the  goal  of  achieving 
room  temperature  MWIR  lasers  and  detectors.  The  relevant 
time  constants  range  from  nanoseconds  to  tens  of  nanosec¬ 
onds,  which  makes  time  resolved  techniques  based  on  me¬ 
chanical  delay  stages  impractical.  This  is  due  to  the  signifi¬ 
cant  variation  of  the  laser  beam  spot  size  and  misalignment 
involved  along  the  delay  path  length.  Alternatively,  while 
fast  MIR  detectors  capable  of  resolving  dynamics  down  to 
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—  30  ns  do  exist,  they  do  not  provide  information  on  the 
radiative  and  nonradiative  processes  involved.  In  this  paper, 
we  report  on  a  novel  technique  to  extract  the  radiative  and 
nonradiative  processes  as  a  function  of  carrier  density  in 
these  lasers  using  nonlinear  correlation  of  photolumines¬ 
cence.  This  was  accomplished  by  using  a  laser  source  that 
can  be  operated  in  both  a  continuous-wave  and  short  pulse 
modes.  Through  fitting  of  the  PL  signal  and  its  nonlinear 
correlation  for  both  laser  modes  of  operation,  the  carrier  life¬ 
time  and  radiative  branching  ratio  were  extracted  as  a  func¬ 
tion  of  density.  The  pulse  repetition  rate  in  the  short  pulse 
mode  acts  as  a  reference  clock  in  the  fitting  procedure  for 
achieving  effective  time  resolution.  We  believe  that  this  is  a 
very  generic  approach  and  can  be  extended  to  various  lumi¬ 
nescing  material  systems. 

II.  EXPERIMENTAL  SETUP  AND  MEASUREMENTS 

The  samples  for  this  study  were  grown  by  solid  source 
molecular  beam  epitaxy.  The  growth  is  initiated  by  a  0.5  /im 
thick  AlSb  layer.  To  allow  the  growth  of  As-free  light- 
emitting  materials,  the  digital  alloy  buffer  layer,  consisting  of 
alternating  thin  (tens  of  angstroms)  layers  of  Alj_vInvSb  and 
Alj_.In.Sb  is  grown  on  the  AlSb  nucleation  layer,  where  y,  z, 
and  thickness  are  chosen  to  provide  the  desired  average  In 
content.  After  the  growth  of  a  suitable  number  of  periods  of 
such  material,  typically  corresponding  to  a  total  thickness  of 
5000  A,  the  thickness  ratio  is  changed  to  correspond  to  an 
increase  in  the  average  In  content,  usually  in  steps  of  about 
9%  or  10%.  There  are  no  threading  dislocations  visible  in 
cross  sectional  transmission  electron  microscopy  images, 
indicating  that  at  least  the  region  close  to  the  top  of  the 
buffer  layer  has  very  few  dislocations.  The  active  region 
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FIG.  1.  80  K  photoluminescence  traces  at  fundamental  frequency  (PLp/) 
and  sum  frequency  (PLp  5um)  spectra  taken  at  70%  of  the  maximum  average 
excitation  power  (Afmax~  1.3  X  1011  cm-2).  The  sum  frequency  correlation 
spectrum  has  a  zero  crossing  at  excess  energy  e0  ~  26  meV  above  the  band 
edge. 

consists  of  three  compressively  strained  (—0.6%) 
In04Ga06Sb/Al02In0  49Sb  wells  grown  on  Al055In045Sb 
buffer. 

We  used  a  Clark-MXR  erbium  doped  fiber  laser  gener¬ 
ating  subpicosecond  pulses  separated  by  7  pp = 24  ns  at 
1.55  /r m  with  —60  mW  average  power.  The  laser  is  addi- 
tively  pulse  mode  locked  using  a  combination  of  polarizers 
and  birefringent  wave  plates.  By  judicious  adjustment  of  one 
waveplate,  the  laser  can  be  reliably  switched  between  mode 
locked  and  cw  operation  of  similar  average  power  P.  Under 
pulsed  excitation,  using  pulses  that  are  orders  of  magnitude 
shorter  than  the  carrier  recombination  time,  the  initial  carrier 
density  n0  is  proportional  to  the  energy  per  pulse  (n0 
ocP-Tpp).  On  the  other  hand,  the  carrier  density  under  cw 
excitation  is  proportional  to  the  average  power  divided  by 
the  density  dependent  recombination  rate  ncw°^PI p(n0), 
where  p(«0)  is  the  inverse  lifetime. 

The  laser  beam  was  split  into  two  equal  beams  that  were 
modulated  using  one  chopper  wheel  at  different  frequencies 
(&>]  =400  Hz,  to2  =  333  Hz).  Both  beams  travel  an  equal  path 
length  to  the  sample  insuring  temporal  coincidence  for 
pulsed  operation.  The  two  beams  were  combined  on  the 
sample  with  a  10  cm  focal  length  lens  producing  a  maximum 
incident  power  density  of  —400  W/cm2.  The  luminescence 
was  then  collected  and  focused  on  an  InSb  cooled  detector. 
The  signal  was  fed  into  two  lock-in  amplifiers  to  simulta¬ 
neously  measure  the  correlation  signal  PL,  sum  at  the  sum 
frequency  tosum  and  PL,- j  at  a  fundamental  frequency  co;  (wj 
or  w2),  where  the  index  i=p  or  cw  for  pulsed  or  cw  excita¬ 
tion,  respectively.  A  long-pass  filter  (>2.4  p,m)  was  used  to 
block  the  stray  laser  excitation  wavelength.  Using  a  Janis 
cryostat,  luminescence  measurements  were  taken  at  room 
temperature  and  80  K.  The  average  laser  excitation  power 
was  changed  over  two  orders  of  magnitude.  Figure  1  shows 
the  80  K  time-integrated  luminescence  spectra  for  PL;)  sunl 
and  PL pj  taken  at  70%  of  the  maximum  excitation  power 


(«  =  0.7/Vmax).  Due  to  the  nonlinearity  of  the  occupation 
Fermi  function,  the  correlation  spectrum  changes  sign  from 
negative  (saturation)  to  positive  (superlinear).  This  zero 
crossing  occurs  at  an  energy  e0  that  is  proportional  to  the 
Fermi  energy. 

We  used  the  zero-crossing  energy  (e0)  in  the  correlation 
spectrum  (PLp  sum)  of  Fig.  1  to  extract  the  excited  carrier 
density.  Following  Chilla  et  al.,  the  time-integrated  spectrum 
at  low  temperature  for  short  pulse  excitation  of  an  initial 
carrier  density  n  and  for  an  effective  density  of  states  g  is 
given  by  PL(«,e) “Tn(n/ge).6  By  solving  for  the  zero¬ 
crossing  energy  of  the  correlation  spectrum  which  can  be 
approximated  as  PLsum(n,e)  =  PL(n,e)-2PL(«/2,e),  we  ob¬ 
tain  e0— 0.25 {nig).  From  Fig.  1,  e0=26  meV  (above  the 
band  edge)  and  using  an  effective  density  of  states  of 
g~  1012  cm-2  eV-1,7  we  obtain  n~  1011  cm-2  and  a  maxi¬ 
mum  carrier  density  used  in  this  work  Amax— 1.5 
X  1011  cm-2. 

The  measured  PL,j  (i=p,  cw;  j  =/,  sum)  versus  carrier 
density  n0  (normalized  to  /Vmax)  are  plotted  in  Figs.  2(a)  and 
2(b)  for  room  temperature  and  80  K,  respectively.  In  this 
case,  the  cw-excited  traces  are  in  effect  plotted  against  a 
scaled  carrier  density  flux  («0/7'pp).  The  solid  lines  represent 
the  curve  fits  used  to  extract  the  recombination  rates  as  will 
be  shown.  A  common  feature  of  the  data  (open  circles)  in 
both  (a)  and  (b)  of  Fig.  2  is  that  PLpy>PLcwy  for  low  exci¬ 
tation  with  PL pf  saturating  at  higher  excitation  becoming 
smaller  than  PLcwj.  For  the  sum  frequency  traces  (solid 
circles),  PLp  sum  turns  negative  at  medium  excitation  levels, 
while  PLcwsum  stays  positive  up  to  high  excitation  levels. 
These  features  result  from  the  difference  between  the  carrier 
density  levels  under  pulsed  and  cw  excitations  and  the  func¬ 
tional  dependence  of  radiative  and  nonradiative  recombina¬ 
tion  rates  on  these  carrier  density  levels,  as  will  be  shown  in 
the  discussion  below. 

III.  MODEL  SETUP 

Luminescence  correlation  is  a  powerful  technique6'8-15 
that  has  so  far  been  underutilized  due  to  difficulties  in  inter¬ 
preting  measurements  involving  time  and  energy  integrated 
correlations  of  temporally  and  spectrally  varying  carrier 
density.6  This  problem  of  correlations  can  be  avoided  by  time 
resolving  the  nonlinearity  of  the  PL  upconversion  signal  in  a 
nonlinear  crystal,  in  a  similar  fashion  to  time  resolving  the 
nonlinearity  of  pump-probe  measurements.16'17  However,  it 
is  desirable  to  avoid  dealing  with  measurements  involving 
several  meters  of  mechanical  delay  for  resolving  nanosecond 
dynamics.  Even  in  the  presence  of  adequate  fast  detectors 
that  would  alleviate  the  problem  of  mechanical  delay  stages, 
the  radiative  and  nonradiative  rates  cannot  be  resolved. 
Therefore  it  would  be  appealing  to  devise  a  simple  technique 
to  extract  these  material  properties  using  time  and  energy 
integrated  PL  signal. 

A.  Quantifying  the  nonlinear  response 

Given  the  sensitivity  of  the  measured  photoluminescence 
(PL)  correlations  to  the  nonlinear  density  dependences,  the 
correct  theoretical  expressions  for  luminescence  have  to  be 
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FIG.  2.  The  measured  time-integrated/energy-integrated  PL, ;  (i=p,  cw; 
y=/,  sum)  vs  the  normalized  carrier  density  nnINmaY  ( A^max  — 1.3 
X  1011  cm-2)  for  300  K  (a)  and  80  K  (b).  The  lower  curves  (solid  circles) 
are  the  sum  frequency  traces  scaled  by  factors  of  10  and  5  for  300  and  80  K, 
respectively.  The  solid  lines  are  the  fitting  results.  The  pulsed  data  and  line 
fits  are  in  gray  while  those  for  cw  excitation  are  in  black. 


PL[«dc  +  njt)]  =  2  («dc)'.  (1) 

;  l- 

By  expanding  the  series  one  finds  that  in  a  lock-in  measure¬ 
ment  the  sum  frequency  term  is  equal  to  PLsum(nmax=n0) 
=  («q/32)o|2PL(«o/2)/ dn2,  and  the  signal  measured  at  one  of 
the  fundamental  frequencies  is  given  by  PLy(nmax=«0) 
=  (n0/4)olPL(«0/2)/ol«.  Therefore  the  measured  lumines¬ 
cence  would  be  proportional  to  PL(n)  only  if  it  happens  to  be 
a  linear  function  which  would  then  imply  zero  luminescence 
correlation.  Similarly,  expressing  the  correlation  term  as 
PL  llm  =  PL( 2n)  -  2 PL (7;)  is  only  correct  for  square  and  cubic 
nonlinearities. 


B.  Luminescence  nonlinearity 

For  short  pulse  excitation,  the  time-  and  spectrally  inte¬ 
grated  PL  is  given  by 

f  rpp 

PL„(«0)  =  c  N(t)p,[N(t)-]dt,  (2) 

do 

where  N(t )  is  the  instantaneous  carrier  density  of  initial  den¬ 
sity  and  n0  and  c  is  a  constant  proportional  to  the  collection 
solid  angle  and  detector  efficiency.  Through  change  of  vari¬ 
ables,  time  can  be  eliminated  from  Eq.  (2)  using  the  rate 
equation  dNI  dt =  -Np{x)  to  obtain 

PLp(«o)  =  ~cf  Pr(N)lp{N)dN.  (3) 

d«o 

Here  we  assume  the  carrier  lifetime  to  be  short  enough  com¬ 
pared  to  7’pp  such  that  there  is  no  constant  accumulated  back¬ 
ground  carrier  density,  which  is  the  case  in  our  sample.  Oth¬ 
erwise  one  has  to  explicitly  take  into  account  the  buildup  of 
a  background  signal  as  a  function  of  the  ratio  between  the 
carrier  lifetime  and  the  pulse  repetition  time.  Fortunately  we 
do  not  need  to  integrate  Eq.  (3)  since  the  measured  quantities 
are  expressed  in  terms  of  derivatives  of  this  integral.  This 
gives 


used  for  reliable  extraction  of  the  inverse  lifetime  p(n)  and 
its  branching  into  its  radiative  pr(n)  and  nonradiative  pnr(«) 
constituents  out  of  the  fitting  procedure.  Even  then,  matters 
are  further  complicated  by  the  confusion  about  what  is 
actually  measured  in  a  lock-in  amplifier  and  what  the  non¬ 
linear  correlation  signal  relates  to.  For  example,  in  a 
modulation  based  technique  a  lock-in  measures  a  quantity 
proportional  to  the  derivative  of  a  function  rather  than  the 
value  of  a  function  as  is  conventionally  assumed.  In  the 
presence  of  two  modulated  pump  beams,  the  excited  carrier 
density  can  be  expressed  as  a  constant  plus  a  time  varying 
sinusoidal  contribution,  ndc  +  nm(t)  =  n0/2+n0l  4{sin(a>lt) 
+  sin(«2?)}>  where  n0  is  the  peak  carrier  density.  In  order  to 
determine  how  the  measured  sum  frequency  signal  relates  to 
the  luminescence  function,  we  expand  the  PL  function  as  a 
Taylor  series  around  the  dc  value, 


PLp/n0)  - 


n0  pr(n0/2) 
4  p(n0/2)  ’ 


(4) 


—  c  22 


1  dPr 
p{n  q/2)  dn 


nJ2 


P,(«q/2) 

p2(«0/2) 


dp 


dn 


V2- 


(5) 


The  cw-excited  luminescence  signal  is  measured  for  the 
same  average  power  used  in  pulsed  excitation.  This  allows 
expressing  the  cw  luminescence  signal  in  terms  of  the  same 
average  carrier  flux  n0ITpp  used  for  the  pulsed  measure¬ 
ments. 


PLcwlVLp)  =  «cw(ho)P/-[«cw(«o)]Lp’ 


(6) 
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«ol 


PLcw,/(V7’pp)  =  c7pP  4 


(  no\ 

dncw 

p, 

IMtJJ 

dn 

tiJ2 


/«o\ 

^Pr(^cw) 

] 

\  2/ 

dn 

V2  J 

(7) 


PLcw,sum(V7’pp)  =  cTPP^\ 


(  no\ 

^cw 

r™U/J 

dn 2 

nJ2 


+  2 


dnr 


dn 


nj  2 


dp,(n  cw) 


dn 


V2 


(8) 


where  ncw,dncwl  dn,  and  cP-ncvll dn1  can  be  obtained  from  the 
relation  ncw(n0) p[ncw(n0)]=n0l  Tpp.  Note  that  the  dependence 
of  Eqs.  (7)  and  (8)  on  Tpp  but  not  in  Eqs.  (4)  and  (5)  acts  as 
a  reference  clock  and  an  anchor  in  the  extraction  of  p{n)  by 
fitting  the  data  to  Eqs.  (4),  (5),  (7),  and  (8). 


IV.  SOLUTION  METHODOLOGY 

The  treatment  so  far  has  been  general  and  applicable  for 
all  luminescing  materials  whether  molecular  or  semiconduct¬ 
ing  and  no  specific  radiative  and  nonradiative  processes  have 
been  assumed.  Solving  the  problem  begins  by  choosing  an 
analytic  expression  for  the  radiative  and  nonradiative  rates 
that  reflects  the  significant  processes  in  the  material  system 
used  whether  Auger,  electron  transfer  quenching,  carrier  dif¬ 
fusion,  etc.  In  this  work,  for  a  semiconducting  material,  the 
decay  rates  were  represented  by 


Pr(«)=PrO+Prl"  +  Pi-2«2> 

(9) 

,  \  Pnrl«+Pm-2«2 

PnrW-Pnr0+  ,  ,  IAT 

1  +  n/Nsm 

(10) 

pin)  =  pr(n)  +  pm[n) . 

(ID 

These  forms  allow  for  simulation  of  linear,  quadratic,  and 
saturating  dependence  on  the  carrier  density  as  needed  for 
realistic  modeling  radiative  and  Auger  processes.  Approach¬ 
ing  degeneracy,  saturation  starts  to  set  in  for  both  radiative 
(pr2  <  0)  (Refs.  18  and  19)  and  Auger  (at  Nsm)  (Refs.  20-22) 
rates.  These  forms  also  allow  for  simulating  the  effect  a  con¬ 
stant  dopant  density  has  on  the  PL  density  dependence.  Some 
of  these  coefficients  could  be  deemed  insignificant  in  a  spe¬ 
cific  temperature  and/or  data  range  as  a  result  of  the  data 
fitting  procedure. 

We  solve  for  the  eight  parameters  of  Eqs.  (9)  and  (10) 
and  the  unknown  parameter  c  of  Eq.  (2)  by  simultaneously 
fitting  the  four  measured  pulsed  and  cw  luminescence  curves 
to  their  corresponding  expressions  in  Eqs.  (4),  (5),  (7),  and 
(8).  In  principle  it  is  possible  to  extract  all  these  parameters 
from  only  one  of  the  four  curves  if  data  could  be  acquired 
with  zero  laser  noise  and  drift  over  a  large  enough  range  and 
with  the  proper  density  of  points.  By  using  all  four  curves  in 
a  simultaneous  fitting,  we  compensate  for  these  experimental 
limitations.  In  order  to  perform  this  simultaneous  fitting,  we 
exploit  the  orthogonality  of  sinusoidal  functions  of  different 


frequencies  to  form  an  orthogonal  superposition  of  the  four 
measured  curves  3L.xp(«)  and  fit  that  to  the  corresponding 
orthogonal  superposition  3{n  ;pr,pnl.)  of  Eqs.  (4),  (5),  (7),  and 
(8)  given  by 


3(n;pr,Pnr)  =  PLP,f+  PLp.sum  sin(10n)  +  PLcw/sin(40n) 

+  PLcw.sum  Sin(  100«) .  (12) 


This  insures  that  the  four  equations,  and  their  correspond¬ 
ing  measured  curves,  contribute  to  the  error  function 
[j(n ;  p, ,  pnr)  -  3exp(n)]  at  different  modulation  frequencies. 
Thus  enabling  us  to  visually  ascertain  whether  the  fitting/ 
minimization  result  is  satisfactory  for  all  four  equations  or 
whether  the  fitting  parameters  need  to  be  perturbed  to  allow 
convergence  to  a  different  set  that  minimizes  the  overall  cost 
function. 

Since  there  are  eight  fitting  parameters  (c,  Nsm,  pri,  pmi, 
/  =  0,1,2),  the  chi-square  as  a  function  of  the  parameter 
space  could  be  a  very  rugged  terrain  with  no  smooth  trajec¬ 
tory  toward  the  absolute  minimum.  Therefore,  it  is  important 
to  guide  the  fitting  program  by  narrowing  the  available  pa¬ 
rameter  space  using  constraints  on  the  forms  of  physically 
acceptable  solutions  to  prevent  the  optimization  routine  from 
being  trapped  in  a  local  minimum  representing  a  subopti¬ 
mum  or  even  a  nonphysical  solution.  One  such  constraint  is 
for  the  total  recombination  rate  to  be  a  monotonically  in¬ 
creasing  function  of  n.  A  more  powerful  constraint  is  formed 
by  forcing  the  optimum  solution  to  satisfy  ratio  of  measured 
PI w  at  two  different  temperatures.  From  Figs.  2(a)  and  2(b) 
the  average  ratio  PLpf  at  80  K  to  PLpf  at  300  K  is  roughly  ~  16. 
Thus  Eq.  (4)  leads  to  the  following  condition  to  be  satisfied 
by  the  optimum  solution  at  both  temperatures, 


PLJbo 

PLJrt 


(13) 


For  experimental  confirmation  we  time  resolved  the  total 
recombination  time  by  measuring  the  correlation  as  a  func¬ 
tion  of  the  time  delay  between  the  two  pump  pulses.  Taking 
into  account  the  nonmonotonic  dependence  (displayed  in 
Fig.  2)  of  the  correlation  on  the  carrier  density,  the  extracted 
carrier  recombination  times  were  5  ±  0.5  ns  at  80  K  and 
0.8  ±0.3  ns  at  300  K.  These  values  are  effectively  an  aver¬ 
age  lifetime  taken  at  maximum  excitation  since  the  lifetime 
starts  shorter  and  increases  with  diminishing  carrier  density. 
The  optimization  routine  was  executed  with  various  starting 
parameters  close  to  and  far  from  the  above  measured  recom¬ 
bination  times.  It  was  found  that  only  one  optimal  solution 
that  is  consistent  with  the  measured  values  produces  the  low¬ 
est  chi-square  fit.  Starting  from  another  value  more  than  a 
factor  of  2  away  from  the  measured  rate  often  results  in  the 
routine  being  trapped  in  a  local  minimum  with  a  chi-square 
value  of  at  least  seven  times  that  of  the  optimal  solution. 


V.  DISCUSSION 

The  extracted  recombination  functions  are  plotted  in 
Figs.  3(a)  and  3(b)  at  room  temperature  and  80  K,  respec¬ 
tively,  along  with  the  experimentally  measured  rates.  The 
recombination  rates  are  extracted  for  carrier  density  ranging 
from  0  to  ndc=n0l  2  since  the  coefficients  of  the  Taylor  ex- 
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FIG.  3.  The  extracted  recombination  functions  are  plotted  for  300  K  (a)  and 
80  K  (b)  as  a  function  of  the  normalized  carrier  density.  The  measured 
recombination  rates  are  indicated  by  the  open  circles. 


pansion  in  Eq.  (1)  are  expressed  in  terms  of  ndc.  Several 
features  responsible  for  the  measured  cw  and  pulsed  nonlin¬ 
ear  responses  are  evident  in  the  solution.  The  radiative  rate 
first  increases  strongly  with  n0  and  then  saturates.  Therefore 
the  ratio  pr{n)l  p{n),  which  is  the  slope  of  PL;,  d,  initially 
increases  and  then  decreases.  On  the  other  hand,  PLCW  ^  is 
determined  by  the  radiative  rate  at  a  lower  carrier  density 
ncw(« o)  rather  than  n0  where  pr(n)lp(n)  increases  strongly. 
That  explains  the  behavior  of  PL,y  as  a  function  of  average 
power.  Similarly,  due  to  the  strong  saturation  of  the  radiative 
rate,  the  PL  develops  a  negative  curvature,  thus  a  negative 
PLp  sum,  at  high  carrier  density  in  the  pulsed  case  but  not  for 
the  cw  case  where  it  is  determined  by  the  radiative  rates  in 
the  low  density  nonsaturating  regime.  The  radiative  branch¬ 
ing  ratio  drops  from  —54%  at  80  K  to  about  3%  at  room 
temperature,  resulting  from  an  order  of  magnitude  increase 
in  the  nonradiative  rate  combined  with  a  reduction  by  a  fac¬ 
tor  of  2  in  the  radiative  rate. 

The  extracted  fitting  coefficients  along  with  their  stan¬ 
dard  deviation  are  listed  in  Table  I  for  both  300  and  80  K. 


TABLE  I.  The  extracted  fitting  coefficients  along  with  their  standard  devia¬ 
tion  for  both  300  and  80  K 


300  K 

80  K 

pro=0.03  ±  0.004  ns-1 

pro = 0.048  ±  0.005  ns-1 

prl  =  (0.08±0.05)/Nmax  cm2 ns-1 

prl  =  (0.5  ±0.08)//Vmax  cm2  ns-1 

pr2=(-0.09±0.05)//Vmax2  cm4 ns-1 

Pr2  =  (-0-56±  0.073) //Vmax2  cm4  ns-1 

PmO=0.93  ±  0.006  ns-1 

pmO=0.097±  0.005  ns-1 

Pnri  =  (21  ±  1.0) /AUx  cm2  ns-1 

Pnri  =  (0.07  ±0.006 )//Vmax  cm2  ns-1 

pm2  is  not  significant 

Pnr2  =  (0.08±0.1)/Amax2  cm4  ns-1 

NSm=(0.4±0.5)Nmax  cm-2 

NSnr  is  not  significant 

Clearly  the  coefficients  are  not  all  equally  constrained.  The 
coefficients  with  the  highest  confidence  are  pr0  and  pm0  for 
both  300  and  80  K,  whereas  those  with  the  least  confidence 
are  Ns  and  pnr2.  Moreover,  the  coefficients  for  the  80  K  are, 
in  general,  more  constrained  than  those  for  room  tempera¬ 
ture.  This  is  reflected  in  the  better  curve  fits  for  the  80  K  data 
shown  in  Fig.  2.  In  addition,  the  standard  deviation  for  a 
particular  coefficient  is  also  influenced  by  the  sensitivity  of 
these  functions  to  slight  perturbations  of  that  coefficient.  Fur¬ 
ther  experimentation  with  different  analytic  forms  for  pr(n) 
and  pm(n)  and  more  efficient  optimization  techniques  will  be 
carried  out  to  maximize  the  sensitivity  to  all  coefficients  in¬ 
volved  and  to  better  optimize  the  solution.  Finally,  care  has 
to  taken  during  data  acquisition  to  guard  against  laser  drift 
and  beam  misalignment  since  the  sum  frequency  signal  is 
very  sensitive  to  beam  spatial  overlap. 

In  conclusion,  we  report  a  novel  approach  to  determine 
the  radiative  and  nonradiative  processes  in  Sb  based  MWIR 
lasers  using  a  nonlinear  correlation  of  PL  with  pulsed  and 
continuous-wave  excited  PL.  Fast  decay  rates  and  their  tem¬ 
perature  variation  were  extracted  using  a  fitting  procedure 
without  resorting  to  ultrafast  detectors  or  to  time  resolved 
techniques  involving  very  long  mechanical  delay  stages.  We 
believe  that  this  is  a  very  generic  approach  and  can  be  ex¬ 
tended  to  various  material  systems. 
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